Photon generation in an electromagnetic cavity with a time-dependent boundary 



CM. Wilson\ T. Duty^'^ M. Sandberg\ F. Persson\ V. Shumeiko\ P. Delsing^ 

^ Microtechnology and Nanoscience, Chalmers University of Technology, S-41296, Goteborg, Sweden, and 
^School of Mathematics and Physics, University of Queensland; St. Lucia, QLD 4072 Australia. 

(Dated: June 15, 2010) 

We report the observation of photon generation in a microwave cavity with a time-dependent bound- 
ary condition. Our system is a microfabricated quarter- wave coplanar waveguide cavity. The elec- 
trical length of the cavity is varied using the tunable inductance of a superconducting quantum 
interference device. It is measured in the quantum regime, where the temperature is significantly 
less than the resonance frequency 5 GHz). When the length is modulated at approximately 
twice the static resonance frequency, spontaneous oscillations of the cavity field are observed. Time- 
resolved measurements of the dynamical state of the cavity show multiple stable states. The behavior 
is well described by theory. Connections to the dynamical Casimir effect are discussed. 

PACS numbers: 85.25.Cp, 42.65.Lm, 05.45.-a 



Photons, in contrast to electrons, do not directly in- 
teract with each other. However, effective interactions 
can be induced when photons interact with a nonlin- 
ear media. In many cases, these effective interactions 
result in so-called parametric processes which are very 
important in widespread technological applications and 
also in fundamental studies of quantum electrodynamics 
[Tl[2]. In the optical regime, the effective nonlinear it ies 
are generally bulk properties of materials or plasmas. In 
the rf and microwave regime, they are often created by 
lumped-element electrical devices [3^ . These types of sys- 
tems have been extensively studied both theoretical and 
experimentally. 

Alternatively, it has been proposed that parametric 
processes can be observed in a system that is essentially 
linear, but where a boundary condition of the electro- 
magnetic field can be changed rapidly in time. For in- 
stance, this is the central theoretical problem in the field 
of the Dynamical Casimir effect (DCE) [4l|5]. One strik- 
ing prediction of the DCE is that real photons can be 
generated out of the vacuum by the changing bound- 
ary condition. The experimental system often imagined 
in describing this vacuum DCE is a cavity with a mov- 
ing mirror. Estimates suggest this is technically a very 
challenging route to pursue, however. Therefore, it has 
been suggested that the observation of photon generation 
starting from a thermal field, which we will call the ther- 
mal DCE, would be much easier, while still demonstrat- 
ing many of the essential features of the vacuum DCE 
[6]. Alternatively, it has recently been proposed that the 
same physics can be studied by instead modulating the 
electrical length of a cavity [7 . 

Regardless of the source of the interactions, the quan- 
tum behavior of parametric systems has been of great 
interest. For example, the quantum dynamics of para- 
metric oscillators (PO) has been described theoretically 
in several different contexts [8l|9]. In fact, the DCE can 
be mapped to the quantum version of the PO in spe- 
cial cases p^. A driven PO can exist in a number of 



stable dynamical states. In a PO with negligible loss, it 
is predicted that quantum tunneling between the states 
is possible. This coherent coupling between the states 
is then predicted to lead to the formation of so-called 
Schrodinger cat states. In POs with more loss, it is pre- 
dicted that quantum tunneling will always be hidden by a 
distinct process know as quantum activation [11 . Quan- 
tum activation is an incoherent process, driven by relax- 
ation (spontaneous emission), that nonetheless leads to 
" over-the-barrier" switching of the dynamical state, even 
at zero temperature. None of these quantum predictions 
have been conclusively observed in experiment. 

To this end, superconducting systems are interesting 
for several reasons: Josephson junctions can be used to 
make the resonators frequency tunable over a wide range; 
various nonlinearities can be easily designed at a desired 
strength; and all of this can be done with very little dis- 
sipation. In addition, superconducting resonators easily 
work at microwave frequencies, implying that operation 
in the quantum regime is possible using standard cryo- 
genic techniques. Parametric devices based on Josephson 
junctions were pioneered in the 80s by Yurke et al. [12], 
but have recently regained interest [TSHIH]. 

In this work, we have studied tunable, high-Q super- 
conducting cavities. They are quarter- wavelength copla- 
nar waveguide cavities fabricated on-chip (Fig.[T]). The 
cavities are made tunable by incorporating a Supercon- 
ducting QUantum Interference Device (SQUID), which 
is used as a parametric inductance. The cavities are 
cooled such that ksT <C huo^ uq being the resonance 
frequency. This means that the average photon number 
in the cavity is well below one, ie., the cavity is in the 
quantum regime. When we strongly drive the magnetic 
flux through the SQUID at approximately twice the reso- 
nance frequency, we observe the generation of microwave 
photons in the cavity (Fig. 

Our cavities are spatially extended systems described 
by a wave equation. The modulated SQUID imposes 
a time-dependent, nonlinear boundary condition on the 



2 





LP = lowpass filter 



FIG. 1: a) A micrograph of a quarter- wavelength cavity. The 
cavity is probed through a small coupling capacitance while 
the other end is terminated to ground through a SQUID. 
Changing the external magnetic field through the SQUID loop 
changes the boundary condition of the cavity, b) Simplified 
block diagram. The cavity is measured using a circulator and 
a cold amplifier. The magnetic field is applied via on-chip 
control lines, one for high frequency and one for DC. 

wave equation. For a static but nonlinear boundary con- 
dition, an exact, analytical solution is already not possi- 
ble, although, approximate solutions can be found in the 
form of coupled oscillators [21] . Treating time-dependent 
boundary conditions is one of the essential aspects of 
the Dynamical Casimir problem. For a general time- 
dependent boundary condition, the system can be re- 
duced to an effective set of coupled modes where both 
the mode frequencies and the couplings are time depen- 
dent [ini [20]. If the cavity is designed such that the 
static mode frequencies are not multiples of each other, 
then we can find an approximate solution when the mod- 
ulation frequency is twice the static resonance frequency 
of one of the modes. In this case, the dynamics of the 
isolated, pumped mode can be reduced to those of a PO 
[10]. Starting from a pure vacuum state in the cavity, 
the appearance of parametric oscillations in this mode is 
then one example of the vacuum DCE. 

Some care must therefore be taken in describing the 
dynamics of the system, taking into account the difficul- 
ties and subtleties mentioned above. That said, we do 
find experimentally that the dynamics of the system can 
be described in a reasonably quantitative fashion by ([T]) 
with the appropriate set of effective parameters. The 



dynamical variable is the canonical fiux of the cavity 
^ = ^oV^/27r = j_^Vdt' where V is the voltage and 
^0 is the flux quantum. 

A PO can be described by the differential equation 0: 

dtt^ + 2Tdt^ + (cj^ + F cos {upt)) Lp^^Lp^ = ^(t), (1) 

where F is the amplitude of the frequency modulation, 
2r = ujq/Q is the linewidth of the resonance and de- 
scribes the damping of the system, and 7 represents the 
dominant nonlinearity of the system, the so-called Duff- 
ing term, ^(t) represents a mean-zero noise term that 
leads to activated switching. 

We derive an expression for F by linearizing the tun- 
ing curve of the cavity cJq = cj^(l + r(x))~^ with re- 
spect to the external flux x = ^^/^o [16 • Here, 
r{x) = Ls{x)/dLi where Lg is the SQUID inductance, d is 
the physical length of the cavity, and Li is its inductance 
per unit length. We find F = 27r tan(7rX(ic)(^o/^&)^^ac 
where the subscripts dc and ac distinguish between the 
static bias point and ac pump amplitude. In our sys- 
tem, 7 arises from the coupling of the current in the 
excited cavity to the SQUID. Following [21], we find 
7 = -uj'^S^/tt with S = -(7r/2)r/(l + r). 

Once the system has been put in the form of ([T]), its 
behavior can be understood in terms of the slow quadra- 
ture variables, qi and (72, where (p{t) = qi{t) cos (ujpt) — 
q2{t) sm{ujpt). In the rotating frame, the dynamics of qi 
and q2 are determined by the metapotential 

9{qi,<l2) = f (9? + 4) + ^^ {ql -<lf) + ^ {if + 4f (2) 

where ft — — cjq) is the normalized detuning, 

C = F l2Vujp is the normalized drive strength, and /3 = 
Z^l^ujp is the normalized nonlinearity. Below the pump 
threshold value of ^ = 1, this potential has only one min- 
imum centered at the origin. Thus the system does not 
oscillate. For small detunings, as the threshold is crossed, 
two symmetric minima develop, yielding two stable, oscil- 
lating states phase shifted by 180°. As the pump is blue 
detuned, we reach a bifurcation point where the two sta- 
ble states merge into a single "quiet" state at the origin. 
If the pump is instead red detuned, another bifurcation 
point is reached where a metastable state develops at 
the origin. The system then has three states, the two 
TT-shifted oscillating states and a quiet state. 

The cavity and SQUID are fabricated in aluminum us- 
ing the standard double-angle, shadow evaporation tech- 
nique. We have previously shown that the frequency of 
such a cavity can be changed much faster than the life- 
time of the photons in the cavity [16] . The samples were 
mounted in a magnetically shielded sample holder in a 
dilution refrigerator with a base temperature of approx- 
imately 20 mK. The measurements of the cavities were 
done using a cryogenic amplifier at 4.2 K that had a nom- 
inal noise temperature of 4 K. The sample was connected 
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FIG. 2: (Color Online) The output of the pumped cavity, a) Two output power spectra referred to the cavity output, for 
different working points. The offset frequency is centered at /p/2. The widths of the peaks are the switching rates between 
different states (see Fig. [s]). b) Integrated output power as a function of fp for Ap = 0.23. The amplifier noise power has been 
subtracted. The predicted linear dependence on detuning from the bifurcation point is clear, c) Integrated output power as a 
function of Ap and fp. d) Theoretical fit from solving equation[l] The color scales are the same in both plots. 



to the amplifier via a circulator mounted on the mixing 
chamber. The signal was further amplified at room tem- 
perature. The resonance frequencies and the Q-values of 
the samples were characterized by coupling a heavily at- 
tenuated probe signal to the cavity via the circulator. For 
the photon generation experiments, no signal was applied 
to the cavity but a flux pumping signal at approximately 
twice the resonance frequency was applied to the fast fiux 
line. The output of the pumped cavity was then recorded 
using a vector digitizer. 

We have studied the pumped cavities as a function of 
both pump amplitude and frequency. In Fig.|2j we show 
the output power of a pumped cavity with a;o/27r = 5.18 
GHz and Q = 8900. The pump amplitude, Ap, is scaled 
such that 20\og{Ap) is the power at the microwave gen- 
erator in dBm. In Fig.|2]d, we show the steady state 
solution to ([l]), namely {(p^) = (v'C^ - 1 - Q)/2p. In 
the range where this solution is imaginary or negative, 
we take the solution {(p'^) =0. To scale the y-axis of 
the theory, we define F = FaA^ and calculate a value of 
Fa = 9.9 X lO^^s"^ using a dc calibration of the pump 
fiux coupling. To scale the z-axis (color scale), we calcu- 
late 7 = 1.9 X lO^^s"^ and estimate the total gain of our 
amplifier chain, which is nominally G = 75 dB. To fit the 
data, we allow Fa and G to vary. We find a best fit for 
Fa = 3.0 X 10^^ and G = 73 dB. We note that since the 
power is proportional to G/7, these parameters are not 
independently constrained. It is not surprising that we 
would find a different value for Fa since the cable loss 
and flux mode structure should be different for the dc 
calibration and the ~ 10 GHz pump flux. Keeping that 
in mind, we conclude that the agreement is good. 

The data is clearly asymmetric with respect to fre- 
quency. This is explained in the following way. At the 
blue detuned bifurcation point Q'^ = +(C^ — l)^^^^ the 



oscillating states vanish and only the quiet state is sta- 
ble. In contrast, at the red detuned = — (C^ — l)^^^^ 
the quiet state emerges, but is only metastable, with a 
occupation probability that is exponentially small com- 
pared to the oscillating states. The oscillations only stop 
when the occupation probability of the quiet state be- 
comes significant. In [9 , it is estimated that this happens 
for Q ^ — 4(^. We therefore introduce by hand a cutoff on 
the red detuned side, although we find that Q = — 3.4^ 
agrees better with the data. 

We can visualize the dynamics of the system by mak- 
ing histograms of the measured quadrature pairs (gi, ^2)- 
The maxima of the histograms then correspond to the 
stable points of the metapotential. In this way, we can 
map out the metapotential ([2|. In Fig.jsj we plot his- 
tograms of (^1,^2) sampled at 1 MHz. In agreement 
with theory, we find that the system has three qualita- 
tively different conditions: 1) one stable state where the 
magnitude is zero, 2) two stable states symmetric about 
zero, 3) three states which combine 1) and 2). We can 
also clearly see the softening of the metapotential just 
before the system bifurcates. 

We observe switching between the different states of 
the system. The switching may be caused by a variety of 
mechanisms including thermal activation, quantum acti- 
vation, or quantum tunneling [Hi |22] . A detailed anal- 
ysis of the observed rates including their dependence on 
various parameters is needed to distinguish these mech- 
anisms. This will be the topic of a future paper. 

We thus find a good agreement between the response of 
a cavity with a time-dependent boundary condition and 
the theory of a PO. This experimentally confirms one of 
the fundamental predictions of the DCE literature [10]. 
Furthermore, the quantitative agreement demonstrates 
that the source of photon generation is in fact the time- 
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FIG. 3: (Color Online) Exploring the metapotential. 2D histograms of measured in-phase, qi, and quadrature, ^2, voltages 
for three different working points, (a) Just below threshold, the metapotential softens in one direction and we observe a noise 
ellipse, (b) Above threshold, the system bifurcates, exhibiting two finite-amplitude oscillating states. These states have equal 
amplitude, but are shifted by 180° in phase, (c) When the pump is red detuned from the cavity, a quiet state can coexist with 
the oscillating states, in agreement with theory. The phase rotations between the three images are instrumental artifacts. 



dependent boundary and not some other mechanism. In 
conclusion, we can interpret this work as an observation 
of the thermal DCE. 

In fact, this experiment has been performed in the 
quantum limit. Thus it is fair to ask whether we can in- 
terpret these results in terms of the vacuum DCE. This is 
essentially a question of whether quantum fluctuations or 
thermal fluctuations initiated the cavity oscillations. We 
cannot conclusively distinguish between the classical and 
the quantum result, since we only measure the steady 
state oscillations which are insensitive to the initial con- 
ditions. In the future, it may be possible to distinguish 
the quantum and classical result by observing the sys- 
tem's transient response. Alternatively, the steady state 
solutions are different if the cavity is removed and the 
SQUID is left to modulate the boundary condition of an 
open transmission line [7 . 
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